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ABSTRACT 

Aiming at comparing different morphological models of galaxy clusters, we use two new methods to make 
a cosmological model-independent test of the distance-duality (DD) relation. The luminosity distances come 
from Union2 compilation of Supernovae Type la. The angular diameter distances are given by two cluster 
models (De Filippis et al. and Bonamente et al.). The advantage of our methods is that it can reduce statistical 
errors. Concerning the morphological hypotheses for cluster models, it is mainly focused on the compari- 
son between elliptical /3-model and spherical jS-model. The spherical /3-model is divided into two groups in 
terms of different reduction methods of angular diameter distances, i.e. conservative spherical /3-model and 
corrected spherical /J-model. Our results show that the DD relation is consistent with the elliptical /3-model at 
lcr confidence level (CL) for both methods, whereas for almost all spherical jfi-model parameterizations, the 
DD relation can only be accommodated at 3cr CL, particularly for the conservative spherical /3-model. In order 
to minimize systematic uncertainties, we also apply the test to the overlap sample, i.e. the same set of clusters 
modeled by both De Filippis et al. and Bonamente et al.. It is found that the DD relation is compatible with 
the elliptically modeled overlap sample at lcr CL, however for most of the parameterizations, the DD relation 
can not be accommodated even at 3<x CL for any of the two spherical /J-models. Therefore it is reasonable that 
the marked triaxial ellipsoidal model is a better geometrical hypothesis describing the structure of the galaxy 
cluster compared with the spherical /3-model if the DD relation is valid in cosmological observations. 

Subject headings: Galaxies: clusters: general — distance scale — X-rays: galaxies: clusters — cosmic back- 
ground radiation — Cosmology: observations — supernovae: general 



1. INTRODUCTION 

Galaxy clusters are crucial to our understanding of 
the universe. For example, galaxy clusters have bee n 
used to derive the Hubble constant dReese et al.l I20QI 
MasonetaJJ 120011: ICunhaetalJ 120071: iFreedman & M adon 



20 10|). to discriminate between different cosmological mod- 
els dMohr et aI] [T99l iThomas et al]fl998t ISuwa et al.ll200l 
iHo et alJl2006l) . to test t he intracluster gas mass distribution 
and temperature profile (ICooravlll998l iPiffaretti et al.l 12003b 
IPuchwein & Bartelmannl2007l:lProkhorov et al.ll201 ll) . and to 
trace out the thermodynami cal history using scaling relations 
amon g cluster observables (Mora ndi et al.ll2007t IShang et all 
2009). These studies often require modeling of cluster mor- 
phology and hence can benefit from a better understanding of 
cluster morphology. 

Hypotheses about an object's three-dimensional proper- 
ties are not easily tested through two-dimensional obser- 
vations (Lee & Suto 120041) . One of the major questions 
about cluster morpholo gy is whether it is spherical or tri- 
axial dFox & Pen 20Q2)- Simulations have shown that dark 
matte r halos are triaxial (Uing & Sutoll2002l iKasun &"E vrard 
2005), and th ere is evidence f rom strong gravitational-lensing 
observations (Gavazzi 2005) as well. Moreover, some ef- 



forts have been made to reconstruct three-dimensional clus- 
ter morphology and correct the projection effect using ob- 
serva t ional data (|Lee & Sutoll2004t IPuchwein & B artelmann 
2006; ISerenoll2007l) . In this paper, we examine the spherical 
and elliptical models of cluster morphology using the cosmic 
distance-duality (DD) relation ( Etherin gtonll 1 9 3 3l also called 
Etherington relation). 

The DD relation plays a fundamental role in observational 
cosmolo gy, covering the a nalyses of galaxy cluster obser- 
vations (Cunha et all 120071) . gravitational-lensing phenom- 
ena (ISchneider et al. 1 ( 1999) and the cosmi c microwave back- 
ground (CMB1 data (iKomatsu et alJl201 lh . The DD relation 
connects different metric distances via 



D L = D A (1 +zY 



(1) 
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where and Da are the luminosity distance and angular di- 
ameter distance respectively, with z being the redshift. The 
DD relation is a general duality in all metric theories of 
gravity, as long as there is no sink or source along the null 
geodesies. 

One way to test the validity of the DD relation is to combine 
the metric distance results from both observations and theoret- 
ical ex pressions in a given cosmological model. lUzan et aTl 
(2004) proposed the idea of testing the DD relation using D\ 
from X-ray surface brightness and S unyaev-ZeFdovich effect 
(SZE, Sunvaev & ZeFdovich|[l972) measurements of galaxy 
clusters. They concluded that there was no significant viola- 
tion of the DD relation fo r a A-C old Dark Matter (ACDM) 
model, de Bernardis et al.l (l2006[) examined the DP re lation 
with Da of 38 clusters from the lBonamente et al. (2006) sam- 
ple. They defined i](z) = D L (z)/ [d a (z)(1 + z) 2 ] and found 
r] = 0.97 ± 0.03 at lcr confidence level (CL). Thus, there 
is no violation for the DD relation in the framework of the 
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ACDM model. With Dl from data compilation of Type la 
supernovae (SNe la), the DP relation was also u sed to con- 
strain cosmic opacity by Avgoustidis et al. (2010). They in- 
troduced a parameter s to analyze the deviation of the DD 
relation in a fiat ACDM model, by assuming that it satisfies 
D L = D A (1 + z) 2+£ , and got t he constraint, e = -0.04+[}^ 
(2cr CD. IHolanda et all d201 ll) offered a method for testing 
the DD relation using WMAP (7-year) results by fixing the 
ACDM model. Particularly, they considered two different ge- 
ometries for galaxy clusters, i.e. elliptical and spherical mod- 
els. Da of these two models were derived by the joint analysis 
of X-ray surface brightness observations plus SZE data. Their 
analysis was based on two parametric representations of tj(z), 
i.e. t](z) = 1 + rj z and t](z) = 1 + t]qz/(1 + Z). They concluded 
that the best-fit value for the elliptical model is close to 770 = 
at lcr CL, whereas for the spherical model, the result is only 
marginally compatible at 3<x CL. 

More robust, another possibility of testing the DD rela- 
tion is via the combination of different sets of observations 
that furnis h both metric distanc es, i.e. Dl and Da, inde- 
pendently. Ho landa et alj (120101) made such kind of cosmo- 
logical model-independent test of the DD relation with Da 
from galaxy clusters observations and Dl from Constitution 
SNe la data. They employed the Da data from two cluster 
models. The first mod el was defined by 25 Da of clusters 
dDe Filippis et al.l [20051) using an elliptical /3-model over the 
redshift interval 0.023 < z < 0.784. The second model con- 
tained 38 Da of clusters (Boname nte et al.l l2006) in the red- 
shift range 0.14 < z < 0.89 observed by Chandra and Owens 
Valley Radio Observatory/Berkeley-Illinois-Maryland Associ- 
ation interferometric arrays. For each cluster, they selected 
one SN la which has the closest redshift to the cluster's, re- 
quiring that the difference in redshift (Az = |z c iuster - zsnI) is 
smaller than 0.005. So a direct test for the DD relation could 
be allowed. They also referred to rj{z) as the aforementioned 
two representations, and obtained 770 = -0.28^ 44 (2<x CL) 
for the elliptical model and 770 = -0.42+°-^ (3cr CL) for the 
spherical model. 

In a subsequent paper, iLi et al.l (1201 ll) also performed a 
cosmological model-independent test of the DD relation, us- 
ing the most recent compilati on which consists of 557 SNe 
la (Union2 compilation, Amanulla h et al.ll2010l). Com pared 
with the Constitution set used by IHolanda et al.l (120101) . the 
values of Az are more centered around the Az = line. Fur- 
thermore, they assumed two more general parameterizations 
for the test of the DD relation, i.e. 77(z) = 770 + r\\Z and 
77(2) = 770 + ?7iz/(l + z). Finally, they obtained the conclusion 
that the DD relation can be accommodated at lcr and 3<x CLs 
for the elliptical and spherical models with 77(z) = 1 + tjqz and 
77(z) = 1 + ?7oz/(l + z), whereas lcr and 2<x CLs by postulating 
two more general parameterizations for the two models. 

There are two aims of this paper: (1) making a cosmo- 
logical model-independent test of the DD relation with two 
new methods using SNe la data from Union2 compilation and 
ga laxy clusters data from two m orphological mod els reported 
by iDe Filippis et al.l (120051) and iBonamente et al.l (120061) . (2) 
testing the intrinsic shape of clusters if the DD relation is 
compatible with present observations. In this paper, two new 
methods to obtain Dl at cluster redshifts are employed: (I) 
fitting the Union2 SNe la data with weighted least-squares 
and interpolating Dl at each cluster's redshift, (II) binning 
SNe la Dl within the redshift range |z c iuster - zsnI < 0.005 
to get Dl at the cluster's redshift. In doing so, we also take 




Figure 1. Quadratic Fit to the Union2 SNe la data. Our best fit and the lcr 
CL are represented by red solid and green dashed lines respectively. The blue 
points with error bars in the upper and lower panel stand for the Union2 data 
and their residuals compared with the fitting function. 



into ac count the asymmetric uncertainties in Bonament e et alj 
(2006)'s Da data. Furthermore, two cluster morphological 
models, i.e. elliptical /3-model and spherical /3-model, are 
tested. The latter model is divided into two groups in terms of 
the different reduction methods of the angular diameter dis- 
tance data, i.e. conservative spherical 73-model and corrected 
spherical /3-model. To avoid limitations of any particular pa- 
rameterizations of the deviation from the DD relation, four 
one-dimensional and four two-dimensional parameterizations 
are applied to the maximum likelihood estimation test. 

This paper is organized as follows. In Section 2 we briefly 
describe the SNe la data and galaxy cluster models. The anal- 
ysis methods and results are presented in Section 3. Section 4 
gives our conclusions and discussions. 

2. SNE IA DATA AND GALAXY CLUSTER MODELS 

In our analysis, the validity of the DD relation is tested 
using Dl and Da results from mutually and cosmologically 
independent measurements. Two sets of Da data from X- 
ray and SZE observations of galaxy clusters are consid- 
ered. In order to get D l, we choose the Union2 sample of 
I Amanu llah et al] (120101) . which contains 557 well-measured 
SNe la. Compared with the "Union" SNe la compilation, 
the new compilation increases the number of well-measured 
distant SNe la by including the events discovered in ground- 
based searches during 2001 and then followed with the Wide 
Field Planetary Camera 2 on the Hubble Space Telescope. 

An important issue is that there are only eight SNe la 
host galaxies for which high-quality Ceph eid distance mea- 
surements are possible (Ries s et al.ll201 ll) . In other words, 
although the sample in Union2 compilation is sufficiently 
large, the calibration of SNe la still suffers from small num- 
ber statistical uncertainty, d ue to lack of SNe la calibrators 
dFreedman & Madore 2010). In order to consider this uncer- 
tain absolute magnitude, we add an uncertainty of 0.05 mag- 
nitudes in quadrature as a covariance among all distance mod- 
uli of 557 SNe la, provid ed by Supernova Cosmology Project 
(lAmanullah et al.ll2010f). 

Both lHolanda et alJ OOlOh and lLiet all (1201 lh employed a 



1 http://supernova.lbl.gov/Union/ 
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moderate redshift criterion, Az = |z c iuster - ^snI < 0.005, and 
selected the nearest SN la for every galaxy cluster to test the 
DD relation statistically. However using merely one lumi- 
nosity distance Z) L from all those available within the same 
redshift range will lead to larger statistical errors J3 Instead 
of using D L of Union2 SNe la directly, two new methods are 
adopted, 

(I) fitting the Union2 SNe la data and interpolating Z) L 
at each cluster's redshift; 

(II) Binning the Union2 SNe la Z) L within the redshift 
range |z c i uster - ZsnI < 0.005 to get D L at the cluster's redshift. 

In realizing method (I), a quadratic fit is applied to all the 
data with errors under the princip le of weighted le ast-squares 
fitting (in other words, x 2 fitting dPress etal.ll 1992b ). The fit- 
ting curve and its lcr CL are plotted along with the original 
data in Figure 1 . 

As suggested by others (Hol anda et al.l 120 lOt iLi et al.l 
12011b . the |Az| < 0.005 criterion is employed to select SNe 
la Z\ data. But in method (II), we take an inverse variance 
weighted average of all the selected data in the following man- 
ner. Assuming that D L , represents /th appropriate SNe la lu- 
minosity distance data (within [Az| < 0.005) with ct Dl . denot- 
ing its reported observational uncertainty, in light of conven- 
tional data reduction techniques by Bevin gton & Robinson! 
(2003, Chap. 4), it is straightforward to obtain 



"Li 



(2) 



A. = 



where Z) L stands for the weighted mean luminosity distance 
at the corresponding galaxy cluster redshift, and erg is its 
uncertainty. 

To get reasonable results for cluster D\, one has to as- 
sume certain cluster morphologies. The structure of clus- 
ters is an essential cosmological probe, as it plays impor- 
tant roles in discriminating different cosmological models 
via the mass de nsity of the universe (Rich stone et al.lfl992t 
ISuwa et aT]l2003h. and constraining the halo evolution mod- 
els dJing & Sutoll2002b . Moreover, different assumptions of 
cluster sha pe affect the measurements of baryon fraction sig- 
nificantly (ICooravlll998l:lAllen et al.ll2004b . Generally speak- 
ing, a robust way to measure cluster D A is through the joint 
analysis of X-ray emission and SZE. The hot intracluster 
medium interacts with CMB photons via inverse Compton 
scattering, causing a change in the apparent brightness of 
CMB and a small distortion of CMB spectrum, i.e. the 
SZE (Sunvaev & Zel'dovich 1972; Birkinshaw 1999)(also see 
Carlstrom et al. 2002J for a recent review). This effect is 
insensitive to the reds hift of galaxy clusters , particularly at 
high redshifts (z > 1) (Carlst rom et al.ll2002b . In calculating 
thermal SZE, on e should consider the relativistic corrections. 
Iltoh et all |l998) proposed a convenient analytical expression 
up to fifth order terms in ksT e /m e c 2 . If the peculiar velocity 
of clusters is non-negligible, namely th e kinematic SZE, th e 
relativistic corrections are presented bv lNozawa et al.l (11998b . 

7 For instance, ther e is a galaxy cluster Abell 1413, at z = 0.142 
iBonamente et al. 2006). Five SNe la in the Union2 SNe la compilation sat- 
isfy the criterion |Az| < 0.005: SN2005fj (z = 0.143), SN2005fw (z = 0.143), 
SN1999bm(z = 0-1441 ), SN2005W (z = 0.145) and SN2005gx (z = 0.146) 
jAmanu llah et al. 2010). Then how to choose one from the 2 SNe la at 
z = 0.143? Using both will improve the statistics. Hence, as a generalization, 
we propose to use all SNe la Dl within |Az| < 0.005. 



Mostly recently, iNozawa et al.l (12006b improve these correc- 
tions to fourth order in ksT e /m e c 2 . 

Two cluster models under different morphological hy- 
potheses are utilized. The first model involves 25 X-ray- 
selected clusters, with m easured SZE temperature decrements 
dDe Filippis et al.l 12005b . A marked triaxial ellipsoidal {3- 
model is reconstructed to describe the cluster structure, as de- 
scribed by 



n e (r) = n e0 



1 + 



2 2 2 2 

e i X l,int + e ^ X 



Z'j.mt + e 3*3,int 



-3/3/2 



(3) 



with x, i nt (i = 1,2,3) defining the intrinsic orthogonal coor- 
dinates aligned with the three corresponding principal axes. 
<?,- (i = 1 , 2, 3) represent the axial ratios (note that one of 
them is unity), and r c denotes the core radius along one of 
the principal axis. As a combination of Chandra, XMM- 
Newton and ROSAT observations, this data set consists of two 
sub-samples, one of which is comprised of 18 galaxy clus- 
ters with 0.14 < z < 0.78 (Ree se et all 12002b . The other 
sub-sample, analyzed by lMason et all (12001b . c ontains seven 
cluster s from the X-ray flux-limited sample of lEbeling et al.l 
( 1996). The second morphological model includes 38 galaxy 
clusters in the redshift r ange from 0.14 to 0.89, provided by 
IBonam ente et aQ (12006b . They modeled the electron density 
profile by a generalized spherical /J-model, 
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where r c i and r C 2 correspond to the two core radii which de- 
scribe the inner and outer portions of density distribution re- 
spectively, and / is a factor between zero and unity. The an- 
gular diameter distance data of every individual galaxy cluster 
with the two models are summarized i n Tabl e 1 . 

As reported by Bona mente et al.l (120061) . almost all an- 
gular diameter distances for the spherical y6-mo del are fol- 
lowed by asymmetric uncertainties. According to lD'A gostini 
(2004), the sources of asymmetric uncertainties include non- 
Gaussianity of the likelihood curve, nonlinear propagation, 
and some systematic effects. And t hus, the real value of p hys- 
ical quantities of interest is biased. Hol anda et al.l (1201 0b and 
ILi et all d201 lb addressed this issue by combining the statis- 
tic al and systematic uncert ainties in quadrature. As stressed 
by IBonamente et all ( 12006b . the modeling uncertainties of the 
angular diameter distances presented in their Table 2, 4 and 
5, contribute to statistical uncertainties. Therefore, it is rea- 
sonable to make t he following corrections and estimations 
dD'Agostinil 12004 equations (15) and (16)), i.e. E(D A ) ~ 
D A +0(A + -A„), cr Dk ~ (A + +A_)/2. We also use the reported 
D A value as the expected value (E(D A )) and the larger flank 
of each two-sided error (max(A + , A_)) as the standard devia- 
tion (cr Z ) A ). Hereafter we name the original spherical D A set 
the conservative spherical /J-model, and the D A set, which has 
been applied a moderate compensation shift to, the corrected 
spherical y6-model. Results of the angular diameter distances 
with conservative and corrected spherical /?-models are also 
listed in Table 1. Separately, we extract an overlap sample 
of clusters that are analyzed by both elliptical and spherical 
models. The result from the overlap sample provides a direct 
comparison of two morphological models of galaxy clusters. 

3. ANALYSIS METHOD AND RESULTS 
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Figure 2. Likelihood distribution functions for elliptical /3-model in four one-dimensional parameterizations. Panel (a) corresponds to the results of SNe la data 
obtained with method (I), and panel (b) corresponds to those with method (II). The inserts of panel (a) and (b) give zoom-in views near the peaks of the likelihood 
functions. 
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Figure 3. Likelihood distribution functions for conservative spherical /3-model and corrected spherical /?-model in four one-dimensional parameterizations. These 
two models differ with how the original data with asymmetric unce rtainties are tr eated. The upper two panels use the original spherical Da set and the lower two 
employ the Da set, which has been shifted in terms of the method of D'Agostini 1 2004). The left two panels correspond to method (I) while the right two method 
(II). 
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Figure 4. Likelihood distribution functions for the overlap sample under elliptical /3-model in four one-dimensional parameterizations. Panel (a) corresponds to 
the results of SNe la data obtained with method (I), and panel (b) corresponds to those with method (II). The inserts of panel (a) and (b) give zoom-in views near 
the peaks of the likelihood functions. 
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It is straightforward to introduce the test parameter 77 ac- 
cording to Eq. (|T), 



D L (z) 



D A {z)(l+zy 



= n- 



(5) 



Then, the 770-associated one-dimensional parameterizations 
for 77 are expressed as follows, 



t](z) = 1 + 770Z, 
?7(z) = l +Voz/(l +z), 
i 1 {z) = l+ m zl(\+zf 
rj{z) = \ -?7oln(l + z). 



(6) 
(7) 
(8) 
(9) 



Note that D A (z) can not be obtained directly from the tech- 
niques of S ZE pl us X-ray surface brightness observations. 
lUzan et all d20M) gave a relation that describes the angular 
diameter distance determined from observations, i.e. Da(z) = 
D^ ustel (z)?7 _2 , which reduces to D A (z) only when T] = 1, i.e. 
the DD relation holds with no violation. Combining this re- 
lation with Eq. (0, it is straightforward to achieve the target 
equation for further statistical analysis, i.e. 



D™(z)(l +zY 



(10) 



where the expression of 77 is checked four times according to 
Eqs. ©-(O. Maximum likelihood estimation is employed to 
determine the most probable values for the parameters, via 

L oc e~* I 2 and 



X 



z 



O7(z)-/7obs(z)) 2 



(11) 



where 77 obs (z) = (1 + z) 



■ rvcluster 



(z)/D L (z) and its uncertainty 



^obs 



cluster 



(z)) 



°"Ol(z) 

D L (z) 



(12) 



Utilizing methods (I) and (II) described in Section 2, we 
are able to determine all corresponding data at each clus- 
ter's redshift. Moreover, conservative spherical y6-model and 
corrected spherical /J-model are retrieved from the original 
data of the r < l OOkpc -cut isothermal yS-model reported by 
iBonamente et alj d2006l) (see Table 1). In Table 1, the number 
of SNe la selected for each galaxy cluster for method (II) anal- 
yses is also listed. For the majority of galaxy clusters, more 
than one SNe la are used. Thus the statistics is improved, es- 
pecially for nearby cluster sample, e.g. Abell 1656. Because 
there is no SN satisfying |Az| < 0.005 from the cluster CL 
J1226.9+3332 (the nearest one is at Az = 0.005), this object 
is excluded, and also removed from method (I) analyses for 
the sake of direct comparison between our two methods. Our 
methods not only avoid double counting of any SNe la data, 
but also take into account all possible clusters data. 

According to the results shown in Figure 2 and Table 2, the 
consistency is clear between the DD relation and the ellipti- 
cal geometry hypothesis, since the DD relation value (/70 = 0) 
is always within lcr CL of the best-fit values of elliptical (3- 
model using both methods. However the best-fit values of 
conservative and corrected spherical /J-models depart from 
770 = at nearly 3cr CL (see Figure 3), especially for conserva- 
tive spherical yS-model (both methods) and corrected spherical 
/3-model (method (II)). 



As to the overlap sample, the difference between the two 
morphological models are more prominent. Figure 4 demon- 
strates good consistency between elliptical /J-model and the 
DD relation, with the best-fit values for all parameterizations 
close to the DD relation value (770 = 0). For the two spherical 
models, the results in Figure 5 show a rather marginal com- 
patibility with the DD relation, which is slightly worse than 
the whole spherical sample results (see Figure 3). Statistical 
results of each sample using the two methods are listed in Ta- 
ble 2 for comparison. 

The one-dimensional parameterizations of the deviation 
from the DD relation, i.e. 77 in Eqs. ©-(UK may be some- 
what restrictive in terms of the degrees of freedom. Here the 
DD relation is tested with two-dimensional parameterizations 
of 77, 



n(z)=n +mz, 

rj(z) = rji +7722/(1 +z), 
l(z) = m +V2Z/(l +zf, 
r](z) = m -?72ln(l +z), 



(13) 
(14) 

(15) 
(16) 



where (771,772) has the value of (1,0) if the DD relation 
holds. It is found that (1,0) falls in the lcr region of all two- 
dimensional parameterizations for elliptical /J-model (Figure 
6), regardless of the method used. However, for the conser- 
vative and corrected spherical /J-models, the DD relation can 
only be accommodated at 3cr CL for both methods, except 
the result (2cr) with method (II) for the corrected spherical 
/?-model. 

The results of the two-dimensional analysis of the overlap 
sample (see Figures 8 and 9) are in good agreement with those 
of its one-dimensional analysis. In fact, the overlap sample 
shows a stronger preference for the elliptical geometry, as the 
DD relation is almost always 3cr away from the best-fit pa- 
rameters for the spherical samples. 

Table 3 presents the minimum values of reduced chi square 
(y 2 . /d.o.f.) for each models with the two methods. Gener- 
ally speaking, it shows the fitting quality of the corresponding 
models to each sample of galaxy clusters. Consistent with the 
results given by the best-fit values of 770 in Table 2, the results 
of ^j^/d.o.f. also favor elliptical model given the validity of 
DD relation, since v 2 . /d.o.f. of this model is much closer to 
unity than those of spherical models. 

In sum, not only does spherical models have larger depar- 
tures of the best-fit values of 770 from the DD relation value 
(770 = 0) as shown in Table 2, spherical models also provide 
poorer fits than elliptical models do as presented in Table 3. 
In consequence, we argue that the galaxy cluster sample mod- 
eled by ellipsoidal morphology is a better understanding of 
cluster's intrinsic nature than those modeled by spherical mor- 
phology. 

4. CONCLUSIONS AND DISCUSSIONS 

In this paper, motivated by the investigation on intrinsic 
structure of galaxy clusters, we have tested the validity of the 
DD relation using luminosity distances from Union2 SNe la 
and angular diameter distances from two cluster morpholog- 
ical models, namely, triaxial ellipsoidal y6-model and spher- 
ical y6-model. In order to obtain more reliable results, two 
sub-samples of spherical yS-model are analyzed according to 
different treatments of the two-sided errors of D A data, i.e. 
conservative and corrected spherical /J-models. Moreover, in 
order to directly compare these two types of morphological 
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Figure 6. Likelihood distribution functions for elliptical /J-model in four two-dimensional parameterizations. As in Figure 2, panel (a) corresponds to method 
(I) while panel (b) method (II). The 1-, 2- and 3-tr CLs are plotted by solid, dashed and dotted lines respectively. The pentagram in each panel stands for the DD 
relation value (1,0), while the big dots in corresponding colors represent the actual best-fit values for four parameterizations. 





Figure 7. Likelihood distribution functions for conservative spherical /?-model and corrected spherical /3-model in four two-dimensional parameterizations. As 
in Figure 3, the upper two panels are from conservative model with lower two corrected model. The left two panels are obtained through method (I) while right 
two method (II). The 1-, 2- and 3-cr CLs are plotted by solid, dashed and dotted lines respectively. The pentagram in each panel stands for the DD relation value 
(1,0) while the big dots in corresponding colors represent the actual best-fit values for four parameterizations. 
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Figure 8. Likelihood distribution functions for the overlap sample under elliptical /J-model in four two-dimensional parameterizations. As in Figure 4, panel (a) 
corresponds to method (I) while panel (b) method (II). Different colors of confidence contours refer to different parameterizations as shown in Figure 6. The 1-, 
2- and 3-tr CLs are plotted by solid, dashed and dotted lines respectively. The pentagram in each panel stands for the DD relation value (1,0), while the big dots 
in corresponding colors represent the actual best-fit values for four parameterizations. 





Figure 9. Likelihood distribution functions for the overlap sample under conservative spherical /3-model and corrected spherical /3-model in four two-dimensional 
parameterizations. As in Figure 5, the upper two panels are from conservative model with lower two corrected model. The left two panels are obtained through 
method (I) while right two method (II). Different colors of confidence contours refer to different parameterizations as shown in Figure 7. The 1-, 2- and 3-<x 
CLs are plotted by solid, dashed and dotted lines respectively. The pentagram in each panel stands for the the DD relation value (1,0), while the big dots in 
corresponding colors represent the actual best-fit values for four parameterizations. The inserts of four panels give zoom-in views near the the DD relation value 
(1,0), which is around the edges of 3-<x CL contours of all likelihood distributions. 
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models as well as minimize systematic uncertainties of the 
test, we also conduct the analysis on the overlap sample, i.e. 
the same set of cluster individuals under both ellipsoidal and 
spherical /3-models for electron density distribution profiles. 

To test the DD relation, it is assumed that t](z) = 
D A lustel "(z)(l +z) 2 /D L (z). In practice, the reduction of statistical 
errors is realized by two methods, (I) fitting the Union2 SNe 
la data with weighted least-squares and interpolating Z\ at 
each cluster's redshift, (II) binning the SNe la data within the 
redshift range |z c iuster - ZsnI < 0.005 to get Dl at the cluster's 
redshift. In methods (I) and (II), the tj{z) parameter is param- 
eterized in four one-dimensional forms, rj(z) — 1 + 770Z, t](z) = 
l+77oz/(l+z), r](z) = l+77 z/(l+z) 2 and?7(z) = 1 -?7oln(l +z). 
Other four more general parameterizations are also consid- 
ered, i.e. 7](z) = J]i + TJ 2 Z, T](z) = 771 + /?2Z/(1 + z). l(z) = 
771 + T l2Z/{i + z) 2 and T](z) = 771 - 772ln(l + z), which are desig- 
nated to describe the t wo-dimensional admissible parameter 
space (N air et al.ll2.Ql ll) . 

Maximum likelihood analysis is used to fit the parameters, 
770 or (771, 772). Our results show that for elliptical /3-model, 
regardless of whether method (I) or (II) is employed, the DD 
relation values (770 = or T)\ = 1, 772 = 0) are always in l<x 
region for all one- and two-dimensional forms. The results 
of elliptical /3-model support the idea that elliptical geometry 
is more consistent with non-violation of the DD relation. In 
the case of conservative spherical /3-model, with both meth- 
ods, it is found that the DD relation value is only marginally 
consistent at 3<x with the best-fit values via one- and two- 
dimensional analyses. The result from corrected spherical /3- 
model is that the DD relation can be accommodated at 3<x CL 
for one-dimensional forms except using method (I), and simi- 
lar results is obtained using both methods for two-dimensional 
parameterizations. In consequence, we can not prove that 
spherical /3-model is compatible with the validity of the DD 
relation. 

Nevertheless, it is noticed that these results may be to some 
extent significantly affected by some individual data points. 
For instance, if cluster CL J1226. 9+3332 is still included in 
method (I) (this object is naturally excluded in method (II) 
for there is no SN la satisfying Az < 0.005), the conclusion 
will be transformed into that the DD relation cannot be ac- 
commodated even at 3<x CL for conservative spherical model. 
This result is tremendously different from the original results 
shown in Figures 3(a) and 7(a). Another example is cluster 
RX J1347.5-1 145. If we exclude this cluster from the over- 
lap sample, and apply our method to this new sample, the 
results show that for conservative spherical model the DD re- 
lation can be accommodated at 2<x CL, and its y 2 . / d.o.f . give 

' A mm ' to 

22.470/15 and 20.204/14 for one- and two-dimensional pa- 
rameterizations respectively. This result is not similar to our 
original results neither. Thus some thorough investigations 
on the influence of every cluster on the global compatibility 
should be paid attention to in future work. 

Furthermore, it is worth making a comparison between the 
two methods. Our results demonstrate that the DD relation 
is compatible with elliptical /3-model at lcr not only for one- 
dimensional parameterizations but also for two-dimensional 
forms, regardless of the methods. There is no discrepancy 
between the two methods for conservative spherical /3-model 
neither; the DD relation can be accommodated at 3<x with the 
two methods. Likewise, corrected spherical /3-model has ap- 
proximately similar outcome, using the fitting and binning 
methods. The results of v 2 /d.o.f. are similar between the 

^ mm ' 

two methods, as presented by Table 3. Accordingly, the two 



methods we used are self-consistent. With our analysis, the 
marked triaxial ellipsoidal model is a more reasonable hy- 
pothesis describing the structure of the galaxy cluster com- 
pared with the spherical hypothesis if the DD relation is valid 
in cosmological observatio ns . 

Boname nte et al.l Q2006) discussed the sources of some sta- 
tistical and systematic uncertainties capable of affecting the 
angular diameter distance measurements. The main origins 
of systematic errors are calibrations of SZE, X-ray absolute 
flux and temperature, while the major statistical uncertainties 
come from SZE observation s , X-ra y spectra and images. In- 
terestingly Bonamente et al. (2006) also pointed out that the 
statistical uncertainty by cluster asphericit y has the largest 
effect , i.e. 15% on Da (see Table 3 of Bonamen te et alj 
2006), but will be av eraged out when a large ensemble of 
clusters is considered (Sulkanen 1999). However, our main 
results tell that in describing the three-dimensional morphol- 
ogy of galaxy clusters, the elliptical /3-model prevails over 
the spherical /3-models. We argue that the current clus- 
ter sample with well measured Da is not sufficient large 
to average out the uncertainty by cluster asphericity. Ac- 
tually, several up-to-date observational work strongly sup- 
port the triaxial ellipsoidal morphology for galaxy clusters, 
e.g. observations in X-ray and str ong gravitatiqnal-l ensing 
dMorandietal.1120101) , SZE imaging (ISavers et al.ll2011l) . etc.. 
Recentlv lNagai et all d2007l) proposed another candidate pro- 
file for the distribution of intracluster medium, known as Na- 
gai model, which is based on the well known NFW dark mat- 
ter profile (iNavarro et al.lll996L 119971) . Although analytically 
non-integrable, this model can also reproduce observed in- 
tracluster medium properties through simulations. However 
recent observation is incapable of distinguishing between Na- 
gai mo del and elliptical /3 -model, due to limited spatial reso- 
lution (Sav ers et al.ll201 ll) . Maybe with more advanced detec- 
tion techniques in the future, people can learn more about the 
intrinsic morphology of galaxy clusters from direct observa- 
tions. 
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Table 1 

Galaxy Cluster Data 



Name 


z 


Mpc 


Mpc 


D A |Cons c 
rtl sphe 

Mpc 


Mpc 


num lovi 


i ei 

num - 

'elhp 


num " 

spne 


Overlap Sample 


MS 1137.5+6625 


0.784 


2479 ± 1023 


5070!]^°, 


5070 ± 1960 


5467.5 ± 1695 


2 


2 


2 


MS 0451.6-0305 


0.550 


1073 ± 238 


1470 -|q 


1470 ± 270 


1500 ± 250 


5 


5 


5 


CL 0016+1609 


0.541^ 


1635 ± 391 


1 220 


1220 ± 220 


1242.5 ± 205 


4 


4 


3 


KA J 134/. 5-1 14j 


0.45 1 


1 loo ± zoz 


1830+^f 


5 1 ± 1 20 


517.5 ±115 


7 


7 


7 


Abell 370 


0.374-f 


1231 ±441 


1830 ±410 


1852.5 ± 395 


2 


2 


2 


MS 1358.4+6245 


0.327 


697 ± 183 




810 ±280 


847.5 ± 255 


3 


3 


3 


Abell 1995 


0.322 


885 ± 207 




1200 ±210 


1237.5 ± 185 


3 


3 


3 


Abell 611 


0.288 


934 ±331 


830! 


830 ± 220 


852.5 ± 205 


5 


5 


3 


Abell 697 


0.282 


1099 ± 308 


770+-'° 
' /u -170 


770 ±210 


800 ± 190 


6 


6 


6 


Abell 1835 


0.252 


946 ± 131 


690+160 

-9Q 


690 ± 160 


742.5 ± 125 


6 


6 


5 


Abell 2261 


0.224 


11 18 ±283 


950 +5Ho 
1560+^°' 


950 ± 300 


980 ± 280 


1 


1 


1 


Abell 773 


0.217^' 


1465 ± 407 


1560 ± 360 


1567.5 ± 355 


12 


12 


12 


Abell 2163 


0.202 


806 ± 163 


73> 

' -220 


730 ± 270 


767.5 ± 245 


5 


5*i 


5 


Abell 1689 


0.183 


604 ± 84 


90Q+1BO 


900 ± 190 


877.5 ± 175 


5 


5 


5 


Abell 665 


0.182 


451 ± 189 


760 +160 

1180 -i 


760 ± 160 


767.5 ± 155 


5 


5 


5 


Abell 2218 


0.1 76-' 


809 ± 263 


1180 ±240 


1195 ±230 


5 


5 


3 


Abell 1413 


0.142 


478 ± 126 


620+™ 


620 ± 190 


642.5 ± 175 


5 


5 


5 



Elliptical /S-model only 



Abell 520 


0.202 


387 ± 141 


5*i 


Abell 2142 


0.091 


335 ± 70 


4 


Abell 478 


0.088 


448 ± 185 


4 


Abell 1651 


0.084 


749 ± 385 


1 


Abell 401 


0.074 


369 ± 62 


4 


Abell 399 


0.072 


165 ± 45 


5 


Abell 2256 


0.058 


242 ± 61 


9 


Abell 1656 


0.023 


103 ± 42 


50 



Spherical /j-model only 



CLJ1226.9+3332 





.890 


MS 1054.5-0321 





.826 


RXJ1716.4+6708 


0. 


.813 


MACS J0744. 8+3927 


0. 


.686 


MACS J0647.7+7015 





.584 


MS 2053.7-0449 





.583 


MACS J2129.4-0741 





.570 


MACS J1423. 8+2404 





.545 


MACS Jl 149.5+2223 





.544 


MACS J131 1.0-0310 





.490 


MACS J2214.9- 1359 





.483 


MACS J2228.5+2036 





.412 


ZW3146 





.291 


Abell 68 





.255 


RXJ2129.7+0005 





.235 


Abell 267 





.230 


Abell 2111 





.229 


Abell 586 





.171 


Abell 1914 





.171 


Abell 2259 





.164 


Abell 2204 





.152 



810+ljj 

2090!] f 

1830! R 
730+200 

3580!p,„ 

1220!|° 

1710™° 

1560+f°° 

1500+™ 

1860!g° 

1990!™ 

760+ 19 ° 

-ISO. 

680+;™ 
560!™ 
1 140+3™ 
720 +3?0 

-280 

740+'™ 

' -220 

460+Tio 



510 ±280 
1580 ± 420 
2090 ± 1070 
1830 ± 430 
730 ± 200 
3580 ± 1620 
1220 ± 340 
1710 ±650 
1560 ±400 
1500 ± 760 
1860 ±420 
1990 ± 470 
760 ± 190 
680 ± 270 
560 ±210 
1140 ±370 
720 ± 350 
740 ± 220 
670 ± 130 
510 ±370 
460 ± 110 



855 ± 250 
1655 ± 370 
2292.5 ± 935 
1845 ±420 
752.5 ± 185 
3865 ± 1430 
1265 ±310 

1770: 
1620: 

1695 : 

1920: 

2012.5 ± 455 

767.5 ± 185 

717.5 : 

597.5 : 
1207.5 : 

772.5 : 

702.5 : 

662.5 : 

570: 

467.5 : 



:610 
:360 
:630 
:380 



:245 
: 185 
:325 
: 315 
: 195 
: 125 
:330 
: 105 



3 
6 
4 

gA 2 

8* 2 
4 

1*3 
1*3 

2 
2 
6 
3 
6 
2 

1*4 
2*5 
2*5 

3 
5 



Not e. — The overlap sample are the galaxy clusters studied by De Filippis et al. 1 2005) using elliptical /?-model and meanwhile by Bonamente et al. 
12006) using spherical ^S-model. 

a Computed estimate for the angular diameter distance reported by De Filippis et al. 1 2005) assuming ellipsoidal geometry for galaxy clusters. 

Calculated angular diameter distance reported by Bonamente et al. 1 2006) under spherical symmetry for galaxy clusters. 
c Conservative angular diameter distance estimate for spherical /?-model calculated in this work as described in section 2. 

Corrected angular diameter distance estimate for spherical /3-model calculated in this work as described in section 2. 
e Number of SNe la selected for each galaxy cluster in method (II) as described in section 2. numlovi, nurnUup and num| sp h c stand for the counts of 
selected SNe la for the overlap sample, the whole elliptical sample and the whole spherical sample, respectively. 

^ The two samples by De Filippis et al. 1 2005) and Bonamente et al. 1 2006) give different redshift values for the same galaxy clusters, and therefore 
the final adopted redshift results are taken from the original observations. The redshift of CL 0016+1609 is from Stocke et al. 1 1991) and those of 
Abell 370, 773 & 2218 from Struble & Rool fl99l . 

g There is no SN la satisfying the criterion, Az — kduxier - 2snI < 0.005, for galaxy cluster CL J 1226.9+3332, and thus it is removed from our analysis. 

In order to avoid double counting of SNe la data as well as consider all possible clusters data in our analysis, each of the five pairs of clusters, 
marked by hi, hj, h$, I14, and /15 respectively, has to be binned within, and the members of each pair have to share the same set of SNe la in method 
(II) as shown. The binning algorithm is inverse variance weighted binning, similar to Eq.[2] For a complete comparison, they are also binned in 
method (I) analysis. 
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Table 2 

Summary of one-dimensional maximum likelihood estimation results: r/o- The rjo below is represented by the best-fit 
value at l-cr Confidence Level for each model using two methods. 



Elliptical /3-model Conservative Spherical /3-model Corrected Spherical /3-model 





Parameterization 


Method (I) 


Method (II) 


Method (I) 


Method (II) 


Method (I) 


Method (II) 


Whole Sample 


V 
V 


= 1 + VOZ 

= 1 + )?0Z/(1 + z) 


-0 063+ 0175 
-0.175 

-0.1 03+ - 242 


-() ()47+ 0178 
-Q178 
-0.083+; 


—0 227 +0 - 102 

-0 339+^'^ 
-0.152 

0.282+»| 2 * 


-0 266+ 0102 
-8-1(8 
-0.396+0- | 

-0 557+<> :22 = 
U - JJ -0,225 
n Q9Q+0.126 


-0 156+ 094 

-a23i+8W 


_0 201 +0094 

-Q.Q94 
-0.297+ 8 ' 42 


V 
1 


= 1 + >?oz/(l + z) 2 
= 1 -?/oln(l + z) 


-0 157+^ : ^ 6 
0.082+H 


-0.132+«;f 
0.064+« : 2 !! 


-0.206 

193+°™ 


-0 412+^ :2 ^ 

248+0H7 
U - ZHO -0.117 


Overlap Sample 


V 
V 
V 
V 


= i + voz 

= 1 + ;;oz/(l + z) 
= 1 + ;;oz/(l + z) 2 
= 1 -;7ofn(l + z) 


00 2i+o.i79 
037+ - 251 

-0.251 

0.064+0- 343 
-0.028+0; 2 ' 4 


()03 j +0.1 82 
o.oji_ 182 

045+° 255 

067+° 348 

-0 037+ - 217 


-0 362 +0 - 145 

UJO -0.14! 

-0 497+0- 2 ® 
" -0.209 

-0.654+;j 
°- 428 tn6 


-0.378+0-5 47 

-0.522+0' 2 !! 

-0 691+ 299 
u.uki _ 0299 

449+ 0178 


-0.287 + °;j| 

n -3S7+0' 94 
-U.J8 /_„ 194 

-0.496+ 0275 

a 337+0*. ?6^i 
U - JJ '-0.163 


-0 306+ 0138 

U - H1 '-0.197 

-0 54 1+ - 278 
""^ -0.278 
361+0 .166 



Table 3 

Summary of one- and two-dimensional maximum likelihood estimation results: minimum of reduced chi square 
C^ 2 ota /d.o.f.). The values below represent the^^/d.o.f. for each model using two methods. 



Elliptical /3-model Conservative Spherical /3-model Corrected Spherical /3-model 





Parameterization 


Method (I) 


Method (II) 


Method (I) 


Method (II) 


Method (I) 


Method (II) 


Whole Sample 


V 


= 1 + voz 


24.745/23 


22.031/23 


68.916/32 


66.265/32 


81.575/32 


77.514/32 


v 


= 1+J7 Z/(1 +z) 


24.696/23 


21.989/23 


68.983/32 


66.422/32 


81.640/32 


77.678/32 


v 


= 1+I70Z/(1+Z) 2 


24.645/23 


21.943/23 


69.370/32 


66.971/32 


81.936/32 


78.151/32 


ii 


= 1 -rj ]n(l +z) 


24.719/23 


22.009/23 


68.903/32 


66.290/32 


81.573/32 


77.553/32 


'i 


= vi + viz 


24.126/22 


21.566/22 


67.809/31 


64.947/31 


80.181/31 


75.829/31 


>i 


= m +»bz/(i + z ) 


24.104/22 


21.576/22 


66.525/31 


63.528/31 


78.866/31 


74.334/31 


V 


= vi + »»z/(i + z) 2 


24.069/22 


21.588/22 


64.748/31 


61.592/31 


76.988/31 


72.241/31 


V 


= ;ji - j/ 2 ln(l + z) 


24.118/22 


21.573/22 


67.150/31 


64.219/31 


79.508/31 


75.066/31 


Overlap Sample 


V 


= 1 + 7/0Z 


13.980/16 


12.725/16 


48.284/16 


46.460/16 


56.947/16 


54.320/16 


v 


= 1 + r/oz/(l +z) 


13.972/16 


12.722/16 


48.837/16 


46.971/16 


57.464/16 


54.800/16 


i 


= 1 + T/oz/(l +z) 2 


13.960/16 


12.717/16 


49.626/16 


47.735/16 


58.170/16 


55.493/16 


>i 


= 1 -%ln(l +Z) 


13.976/16 


12.724/16 


48.537/16 


46.689/16 


57.186/16 


54.536/16 


V 


= VI + V2Z 


13.832/15 


12.683/15 


42.376/15 


41.302/15 


50.079/15 


48.419/15 


v 


= VI + V2Z/(1 +z) 


13.811/15 


12.673/15 


40.559/15 


39.517/15 


48.033/15 


46.418/15 


V 


= VI + V2Z/(1 + z) 2 


13.803/15 


12.664/15 


38.222/15 


37.195/15 


45.378/15 


43.785/15 




= Vi - V2ln(l + z) 


13.819/15 


12.678/15 


41.433/15 


40.378/15 


49.018/15 


47.384/15 



